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Abstract 

We consider a family of four-dimensional non-linear sigma models based on 
an 0(5) symmetric group, whose fields take their values on the 4-sphere S . 
An SO(4)-subgroup of the model is gauged. The solutions of the model are 
characterised by two distinct topological charges, the Chern-Pontryagin charge 
of the gauge field and the degree of the map, i.e. the winding number, of the S 4 
field. The one dimensional equations arising from the variation of the action 
density subjected to spherical symmetry are integrated numerically. Several 
properties of the solutions thus constructed are pointed out. The only solution 
with unit Chern-Pontryagin charge are the usual BPST instantons with zero 
S 4 winding number, while solutions with unit S 4 winding number have zero 
Chern-Pontryagin charge. 
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1 Introduction 



The model considered in this work is described by the Lagrangian on R 4 

+ ^\D,r x D u <f) b x D p <p c \ 2 + V(<j) 5 , cosu) , (1) 
in terms of the S 4 valued fields (fi a = (0",0 5 ), (a=l,2,3,4) satisfying the constraint 

= i , 

and the SO (A) gauge connection Ajf^ with curvature Fjfifl. The covariant derivatives 
in ([I]) are defined by 

D,r = v + 4 Q V, d^ = 8^. (2) 

(The brackets [..] imply antisymmetrisation of indices throughout.) 

The Lagrangian ([[]) differs from that of the various models considered in [JTj] , in the 
physically important^] respect that the kinetic term quadratic in the S A field <ft a , which 
was absent there JI|, is present here. As a result of the Derrick scaling requirement, 
the system (P features also a kinetic term sextic in the S 4 valued field. (|TJ) differs 
from that in [II] also in the presence of the generic potential V((j) 5 , cos a;) 

V((p 5 ,uj) = A(coscj - <f) n , < uj < ix , n = integer , (3) 

whose role it is to fix the asymptotic value of the field a , rather like the pion-mass 
potential in the usual 3 dimensional 0(4) Skyrme model. Like in that case |2[], this 
term serves only the purpose of fixing the asymptotics, and will be considered only 
in this limited context. Moreover, anticipating our conclusions in Subsection [2.1| , 
namely that only for uj = is it possible to construct finite action solutions, the 
relevant V((f) 5 ,uj = 0) is 

y(0 5 ,o; = O) = A(l-0 5 ) , (4) 

in analogy with the pion-mass potential. Notwithstanding, in Subsection |2.2| 
we have described a model differing from ([!]) and characterised by nonzero u (in 
particular u> = |), which can support topologically stable finite action solutions. As 
will be seen there, the Lagrangian of such a model does not feature the usual YM 
term and is therefore not studied further here. 

Our topologically stable finite action solutions are interpreted as instantons, al- 
though the latter are not always characterised by the second Chern-Pontryagin den- 
sity, but also by the S A — > S 4 degree of the S 4 valued field (which may or may not be 
integer) . 

1 Without a quadratic kinetic term, it is not possible to infer from finite action conditions, that 
that the matter field (in this case the S 4 valued field) becomes asymptotically a constant and hence 
consistent with vacuum field, ft would then be impossible to interpret the resulting topologically 
stable finite action solution as an instanton. 
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By adapting the methods formulated in || used in establishing the topological 
lower bounds given by the degree of the map S 4 — > S 4 , a suitable such lower bound 
can be established for the action ([]]), for all positive values of Ai, A2, A3, and Ao- This 
will be established in Section Q. In Section |, the system ([[]) will be subjected to 
spherical symmetry, and thenceforth we will restrict to the study of the ensuing one 
dimensional equations numerically. The results of our numerical work is presented in 
Section f|, and Section | is devoted to the summary and discussion of our results. 

In the work of |1|], which employs the above model ([[]) with Ai = A3 = 0, it 
was found that the action of the main 5*0(4) gauged 0(5) model decreased with the 
Skyrme coupling A2 and exhibited a bifurcation at a value very close to and above 
the action of the (pure) Yang-Mills (YM) instanton Q. Thus beyond a critical value 
of the coupling A 2 , the system did not support a finite action solution, and more 
importantly the action could be made smaller by decreasing A 2 . Our main aim in the 
present work is the verification of these two properties of the solutions when Ai > 
and A 3 > 0. We have found that these properties of the solution persist, namely that 
instantons can be constructed for values of A2 (holding A3 constant) up to some value 
AJf, and for values of A 3 (holding A 2 constant) up to some value Ajf. Also, the actions 
of these instantons decrease with A2 and A3 respectively consistently with the Derrick 
scaling requirement. We find surprisingly that the action at A2 = is nonzero, inspite 
of the vanishing of the topological lower bound at that point. We defer discussion of 
the possible physical significance of these properties to Section ||. 

Another objective here is to probe the nature of the topological lower bounds. In 

we were exclusively concerned with lower bounds stated in terms of the degree 
of the map of the S 4 valued field, such that the corresponding solutions supported 
vanishing Chern-Pontryagin charge. Here we attempt to construct numerically, in- 
stantons with nonzero Pontryagin charge when the S 4 field is nontrivial. We do not 
find such solutions and offer an analytic argument to support their nonexistence. 

The boundary conditions employed for the S 4 valued field are 

lime/) 5 = —1 , lim 6 5 = cosu , (5) 

but we will be restricting our numerical investigations to the u = case. For the 
gauge field we will adopt the vacuum behaviours 

lim = , lim = qg d^ 1 . (6) 

q = leads to zero Pontryagin charge. Pure-gauge q = 1, and half-pure-gauge q = \, 
both lead to nonzero Pontryagin charges. The half-pure-gauge case q = \ pertains 
to uj 7^ 0, which we do not study numerically. The pure-gauge cases q = and q = 1 
both pertain to uo = 0, and are studied numerically. It turns out that instantons with 
nontrivial (S 4 valued) matter field can only be constructed for the q = case. 
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2 Lower bounds 



The work in this Section follows very closely that in || and that in Section ^ of |T|]. 
The analysis in |3], [l| was adapted to the definition of topological charges presenting 
lower bounds on the actions, of systems supporting asymptotics with u = in (Q). 
Here, we extend this analysis to include values of < uo < n. In the generic case 
therefore, we would expect solutions with a fractional analogue of the Baryon number 
in 3 dimensions 0, while in the limiting case || [l] u — this will be the degree of 
the map of the S 5 valued field taking an integer value. 

We start with the definition of the winding number density 

Qo = ^e^ pa e abcde d^dA h d p <p c d a <p d <p e , (7) 

which is inadequate for our purposes here since it is gauge variant, and its gauge 
invariant version 

Qg = ^e^e^D^D^D^D^ , (8) 

whose volume integral cannot be evaluated by stating the asymptotic conditions, i.e. 
it is not useful topological charge density. 

The volume integral of (0) can indeed be evaluated by stating the asymptotic 
conditions (|5|), such that for uj = the value of this integral is integer, and is fractional 
for w ^ 0. The normalisation ensures that in the spherically symmetric case this is 
the unit charge, or winding number. The task is to define a suitable density which 
is (a) gauge invariant, and (b) its volume integral equals the volume integral of the 
density g , (0). 

To this end we find the relation between the two densities (|7|) and (§), in suitable 
form, such that all gauge-variant terms appear as total divergences. Thus, 



Qg = Qo + df, U 5 dA p + fi M 
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a(3-yS 



VD^D^F^ - \^ (l - i(0) 2 ) F$ Ff a 



(9) 



where 



= -^^^e a ^ 5 A a r (2d p <^ + {A p( f>y) (10) 

= iJ^w^ 75 (l " \(<P?) [d p Af - \{A p A a r^ } . (11) 

The volume integral of the (gauge- variant) total divergence terms in (||) can, 
after conversion to two surface integrals, be evaluated using the asymptotic values 
d^) and (H). If these surface integrals vanished, then the volume integral of the 
remaining gauge invariant terms could be expressed in terms of the volume integral 
of £>0; namely the (integral or fractional) winding number, leading to the definition of 
a gauge-invariant topological charge density. 
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This can be checked by substituting the spherically symmetric Ansatz ( |3T]) and 
(|32|) in ( p~0|) and (|TT|). One then sees immediately that the density (|T0D vanishes 
asymptotically by symmetry, while the density ([11]) yields a non-vanishing contribu- 
tion to the corresponding surface integral, which depends on the asymptotic param- 
eter to. Thus the definition of the charge density to be given below is tu-dependent. 

Making use of the relation 

2 



■ a/37<5 pa/3 p7<5 _ a/37<5a \ A-fS 

I ^pvpa^ L pv L pa c pvpa& ^p 1 ^>< 



a/3 



:(ApA a ) 



a3 



(12) 



one can add and subtract the gauge-invariant density 

1 , 



cosu; coscj 



cos u e 



- pvpcr'- 



a/S'yS jpa/3 jjryS 



pv pa 



to @, such that the non- vanishing surface contribution of f2 M is cancelled. After some 
rearrangement, the natural definition for the charge density q pertaining to a system 
characterised by the asymptotic parameter uj is 



Q= Qg 



64vr 



2 C W 



e ^lS\Ma D(t) P D( Lh F1 S 



pa 



f 1 



i 



OOSa.' i fOSa' — - COS 2 Co" 



pa/3 p7<5 I 
pv pa ( ' 



(13) 



which is the manifestly gauge-invariant definition that is employed in establishing 
(Bogomol'nyi like) lower bounds, and which is equivalent to the definition 



Q = Qo + dp (0 5 9„n„„ + 
in which Vt u ^ is defined by ( |T0| ) while fi M is 



a. 



128vr 2 ^ 



0° 1 



coscj cosu; 



cos 2 



(14) 



x 



(15) 



For the spherically symmetric fields (|3ll)-(l32|), x^Cl^ vanishes asymptotically, and since 
we already know that Q Ufl vanishes asymptotically, we see that the volume integral of 
( |T4j ) equals the (generic fractional) winding number. The (topological charge) density 
is gauge-invariant, and its volume integral is just the winding number of the S 4 valued 
field. 

Not surprisingly the definition (p~3|) for the two most prominent cases g(u = 0) 
and g(u = ~ ), simplifies somewhat. After some manipulations one has 



g(0) 



Qg + 



Qg + 



64tt 



2 fiupa 



64tt 



2 u pvpa L 



r a/?7<5 j 5 



pa 



(16) 
(17) 
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In the next two Subsections 27T and [2.2| we shall analyse models whose actions are 
bounded from below by the topological charges corresponding to (|16|), pertaining to 
uj = 0, and (|T7|), pertaining to u — ~. 



2.1 Model with u = 

The relevant gauge-invariant charge density in this case is (|T6|) . The first term, q g 
can be reproduced as a consequence of the inequality 



>o, 



leading to 



^l^ri 2 + ^|D[^ [p c A,]]0T > k x kIq g 



d\2 



where k± and k 3 are constants with the dimension of length. The next term in (|TB|) . 



is reproduced by the inequality 



> 



expanding which and adding the appropriate positive terms to the left hand side 
yields 

+ 4\D[^ a D,]<p b \ 2 > nlnle^e^^D^D^F^ . (19) 
Finally the last term in (|TE|) can be reproduced by adding the two inequalities 



5\2 



k 5 pr a P fry 6 

ixv 1 2J2 P°" 



> 



> 



and then adding suitable positive quantities to the right hand side to yield 

1 



2j ' v 4<~ij,vp<j 



a/37<5 [(0 5 ^ 2 



1 1 pa/3 p-yS 
J uu pa ' 



(20) 



and (|20|) , all have the dimension of 



The constants k , in m, K3, and in (|18|) , ([ 
length. 

Adding (0), (|19D and (|C])> we en d up with an inequality whose left hand side, up 
to some redefinitions, is precisely the system ([I]), without the potential term V. This 
Lagrangian is bounded from below by the right hand side, which will be a topological 
bound if the latter coincides with the topological charge density ([ll]) (up to a constant 
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multiple). This turns out to be the case, provided that the constants /to, in K\, k%, 
and k satisfy the following constraints 

with the constant multiplying g (|l6|) , equal to kik^. Thus the action (before redefining 
the constants) is bounded from below as 

S > N , (22) 

where N is the winding number. The action S is the four-volume integral of the 
Lagrange density 

£ = (4 + 4)\F [ ^ ] \ 2 + 4\D,<f> a \ 2 

+4\D^ a x D u <f) b \ 2 + 4\D^ a x D u <p h x /J p c | 2 , (23) 

subject to the restrictions which is up to some redefinitions coincides with (]I|) 
without the potential term V. 

It is seen from (p22"l) that the condition that this lower bound be nontrivial is that 
neither one of K\ and K3 should vanish. It could be thought that this means K2 can be 
set equal to zero without violating this bound, but from the first member of the con- 
straint equations ( pl|) we see that this is impossible. We conclude therefore that the 
lower bound remains valid only as long as all the constants (k , k±, K2, ), 
and hence also the couplings (Ao, Ai, A2, A3, ), remain positive and nonzero. 

As will be shown in Section |3|, the constants Ao and Ai can be scaled away leaving 
only two independent coupling constants A2 and A3, both of which have to be positive 
and greater than zero if the lower bound (|22|) is to be preserved. 

Before proceeding to the next Subsection, we note that the Lagrangian ([!]) is 
not unique in being bounded from below by the topological charge density ([161) . 
Rather, it is the simplest system motivated by the requirements that it features the 
YM term and the quadratic kinetic term of the scalar field. An inspection of the 
spherically symmetric restriction of the YM term, eqn (|33]) below, implies the finite 
action condition on the gauge field function a(r) 

lim a(r) = ±1 , 

which in the language of the asymptotic conditions (||) means that q = and q = 1, 
respectively. In the first case, the Pontryagin charge vanishes, while in the second 
case it is equal to 1. 

2.2 Model with u = f 

The relevant gauge-invariant charge density in this case is (PD. Unlike in the pre- 
vious Subsection however, here we do not proceed straightforwardly to construct the 
simplest density which is bounded from below by (|T7|). 
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The reason is that when w ^ (as in the case case at hand with u — |) the gauge 
group 50(4) breaks down to 5*0(3) at infinity. This can easily be seen by rotating the 
asymptotic field (f> a (a = 1,2,3,4) of length |0 a (oo)| = sintu, to the constant vector 
field along the X4 axis, by means of an appropriate 50(4) gauge transformation. The 
effect of this transformation on the so(4) gauge connection is, that it develops a line 
singularity along the 24 axis, and its non- vanishing components then take their values 
in the residual so(3). We do not give the details of the passage to this 'Dirac gauge' 
here, because this has been given in detail previously in Refs. || [7| , in the context of 
the 50(4) x U{1) gauged Higgs modelf] 

The relevant information that follows from the preceding discussion is, that the 
asymptotic connection field A" 6 = (A^ b , Af 1 ), (a = a, 5) in the Dirac gauge decays 
exactly as -, and its only non- vanishing component is 

r(l — £ 4 ) J 

which takes its values in so(3). 

It follows that the corresponding asymptotic curvature field has the Coulomb 
decay and as a consequence the integral of the YM action density will diverge 
logarithmically in four dimensions. This simple fact can also be seen by inspecting the 
spherically symmetric YM action density in (|33|) . Thus, in constructing the density 
bounded from below by the topological charge density QlTD, it is not legitimate to 
employ the usual YM action density. 

The remedy is to use instead the YM density constructed from the antysym- 
metrised product of two curvature two-forms, namely 

I Tpabcd 1 2 _ I / j?a\p j-,7<5] . plS-y P /3]ax |2 
I^Vpal - + b iAv b p<j\ )\ ■ 

This term arises naturally in reproducing the last term in the charge density fllTD - 
To reproduce the second term in the charge density, it is not legitimate to make 
use of inequality ( |T9D since the latter features the usual YM density. Given that for 
the instanton (vacuum) interpretation of the solution we need to have the quadratic 
kinetic term of the scalar field, this necessitates the appearance of the term 

KPaW ■ 

Finally, to reproduce the first term, g@, in (|i~7|), the most economical option is to adopt 
the inequality ([18]). (This avoids the introduction of the additional and unnecessary 
term {D^^D^^l 2 in the Lagrangian.) 

Following the above arguments, we write down the three topological inequali- 
ties corresponding to (fTcf) , (|19|), (|20D for the present case with uj — |. With (|l"8l) 

2 The analysis here is identical to that in |g with the 4 component field <fr a here replacing the 
Higgs field $ = 757^^ of || 0|. Indeed this is the case for all d-dimensional (d > 3)SO(d) Higgs 
models with d- vector Higgs fields [||, of which the most familiar is the Wu-Yang monopole in d = 3. 
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(24) 



unchanged, we just give the second two 



Rl\D [fi r\ 2 + 4\F^Dx]ff > ^iRle^e^^D^D^F^ (25) 
-4\F^t\ 2 + r\^f > TR^e^e^^fF^F^ , (26) 

where the constants R±, R3 and R4 all have the dimension of length, while r is dimen- 
sionless. 

Adding (|T8|), (p5[) and (|26|) results in an inequality whose right hand side can be 
identified (up to a numerical factor) with the topological charge density (|TTD, provided 
that 

RiR\ = 3kikI , 6TK4 = . (27) 
The resulting topological inequality bounding the action from below, analogous to 

S > N , (28) 
in which the action S is the four-volume integral of 

abed 1 2 1 —61 77>a6 7-> ±c\2 , 6| n J,a . . 7~> jfe . . r-> J.c|2 



^ = ^ICp"I + R b 3 \F^D x] cf>r + 4\D^ a >< x D i 



pi 



+ (kI + kDID^ + t 2 ^) 6 , (29) 

subject to the constraints (p7|). Note that the potential @ with = and n = 6 
appears quite naturally in (^9|), and in this case its presence is mandatory if the lower 
bound on the action is to be preserved. 

Because (p9|) does not feature the usual YM term besides the F A term, it is not 
likely to be of any physical interest. Hence we do not analyse it numerically. We note 
that in the case of the 50(4) x U(l) Higgs model || 0, which also features the F A 
YM term to the exclusion of the usual F 2 , it could be argued that at high energies 
that system reduced to a conventional 5*0(4) x U(l) Higgs system 

Ti(\ 2 F 2 u + X l D ll <t> 2 + X ^ 2 + 7 ] 2 ) 2 ) , (30) 

where the constant rj is the VEV of the Higgs field, and all fields are antihermitian. 
In other words, the 50(4) x U(l) Higgs model was interpreted as the low energy 
effective action of (0). Unfortunately, we do not have such an interpretation for the 



system (|29f) . 



3 Spherical symmetry 

The spherically symmetric Ansatz employed is 

A^ = ^l(S^-8^) (31) 
(f) a = sin f(r)x a , 5 = cos f(r) (32) 
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As explained in Section [2.2| , we will restrict our numerical analysis to the case of 
uj = 0, and hence give the spherically symmetric reduction only of the terms in the 
system ([[]), or pertaining to u — 0. 

Substituting fl31~|)-(|3~2"D into the component terms of ([D we have, for each term 



7>a/3|2 



\D^ a \ 2 



12 



- + 



a 2 -r 



./ 



/2 



'a 2 sin 2 /" 



Lb\2 



( a 2 sin 2 f \ r „.n fa 2 sin 2 f ' 



\D^ a x D„ 



x D^l 2 



36 



'a 2 sin 2 /' 



3/ 



'a 2 sin 2 /" 



(33) 

(34) 
(35) 
(36) 



In the following we will study the classical solutions of the model ([I]) and characterize 
them by the classical action S defined by means of 



S 



111' 



d xC 



(37) 



The reduced 1-dimensional Lagrangian is r 3 times the sum, with the appropriate 
numerical coefficients in ([!]), of all the above 4 terms. We do not display this one- 
dimensional Lagrangian, nor the ordinary differential equations that follow. 
The asymptotic values of the function /(r) corresponding to (||) translate to 

lim/(r) = Ti , lim /(r) = u , (38) 



with uj = 0, while the asymptotic values of the function a(r) for the cases q = and 
q = 1 translate respectively, to 



lim a(r) 
lim a(r) 



lim air) 

r— >oo 

lim afrl 



-1 . 



(39) 
(40) 



Let us first point out that the embedded charge-one-BPST-instanton solutions 

k 2 — r 2 



aBPST{r) 



k 2 + r 2 



f(r) = nn (everywhere) 



(41) 



(k is a real constant, n is an integer) exist irrespectively of the values of and 
leads to Sbpst = f , corresponding to the action of the charge-one- instanton solution 
of the Yang-Mills theory 0]. Here we are interested in classical solutions with non 
constant /(r). 

The number of four coupling constants can be reduced to two by using the follow- 
ing scaling argument. Transforming r — > err, we have 



S(Xi, A 2 , A 3 , A ) = ^(Aicr 2 , A 2 , A 3 <r 2 , A ) = \i<j 2 S 1, 



A 3 



A, 



A f 



Axa 2 ' A lC r 4 ' A lC x 2 



(42) 
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Choosing cr 2 = ^ this gives 

S(\i, A 2 , A 3 , A ) = XoS (l, ^, . (43) 

In the following we will make use of the above scaling property and set Ai = Aq = 1. 



4 Numerical results 

We have studied numerically the solutions of the classical equations associated with 
(HI) for the uj = model, restricting to the one dimensional spherically symmetric 
fields given by (|3~3|)-(|36|). Most of the work is carried out with the potential (|4]) 
decoupled i.e. with A = 0. 

In [|l| the above equations have been studied in detail in the case Ai = A3 = 0. 
Here we want to study the classical solutions for non-vanishing Ai, A2, A3. Using the 
standard Derrick scaling argument, it is easily seen that regular classical solutions will 
exist only if the coupling constants Ai, A3 are both nonzero. On the other hand, the 
topological lower bound ( p2|) derived in the previous Section states that in addition 
to Ai and Ao (which we have already set to Ao = Ai = 1 by scaling), both A2 and 
A3 must be positive and nonzero. On the basis of the last statement, there is no 
guarantee that the solution persists when A 2 vanishes, even though this is consistent 
with the Derrick scaling requirement. 

As a result of our numerical studies, we have learnt that with the asymptotics 
(|39|) , the solution persists when A2 vanishes. In this case there remains only one 
coupling constant to vary, A3, which is a simpler case to study and this is presented 



in Subsection [Q| . In Subsection |4T2| , again with the asymptotics (0), we study the 
cases where A2 is varied for fixed nonzero value of A3, and also where A3 is varied 
for fixed nonzero value of A2. These families of solutions all have Pontryagin charge 
equal to zero. In Subsection ( |4.3| ), we present the results of our numerical search for 
solutions with unit Pontryagin charge and nontrivial scalar field, with asymptotics 
flioD . The result is negative, and we have supplied an analytic construction in support 
of the nonexistence of such instantons. 



4.1 Solutions with a(0) = a(oo) = 1 and A2 = 

With these boundary conditions, the Chern-Pontryagin charge would be zero and the 
topological lower bound would be stated in terms of the degree of the map only. 

Integrating the equations for small values of A3 we were able to construct solutions 
with 

a(0) = l , a(oo) = l , /(0) = tt , /(ao) = (44) 

The profiles of the functions a, / of this solution are presented in Fig.l for A3 = 0.425 
by the solid lines. In the limit A3 = the classical action tends to zero and the function 
a{r) tends to a constant : a(r) = 1. When A3 increases, the function a{r) develops a 
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local minimum (at r = r m ) which becomes progressively deeper as indicated in Fig. 
2. The general dependence of r m on A3 is also reported on Fig. 2. At the same time 
the classical action of the solution increases with A3, and this is illustrated in Fig. 3. 




This situation persists up to a critical value of A3, say A3, and numerically we 
found 

A^ « 0.42661 . (45) 

Corresponding to this critical value, we find r m « 0.041, a(r m ) « 0.034 and S ~ 
1.3345. In particular the value of the action is slightly higher than the value 4/3 
corresponding to the action of the instanton solution of Yang-Mills theory M. 

In fact, a large part of this action comes from the Yang-Mills part of the Lagrangian 
and the contribution due to the S 4 valued matter field is rather tiny (less than one 
percent) because, as indicated by Fig. 1, the function f{r) becomes very steep in the 
region where the function a(r) has its minimum. 

We found no solutions for A 3 > A3; however, a careful analysis of the equations 
strongly suggests that a second branch of solution exists, as illustrated on Fig. 2. As 
far as the classical action of the two branches is concerned they terminate into a cusp 
at A3 = A3, in a way very similar to Fig. 9 of |J. 
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Figure 2 



0.12 




0.44 



X 3 



Figure 2: A3 dependence of r m and a(r m ) for the two branches near the critical value 
of A 3 when A 2 = 0.1. The global dependence on r m is displayed in the window. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 3: A 3 dependence of the action for A 2 = (solid line) and A 2 = 0.1 (dotted 
line) . 
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As suggested by Fig. 2, it is very likely that when A 3 decreases to below A3 
on the second branch, the minimum a(r m ) has a tendency to approach zero while 
the derivative f'(r m ) of the function f(r) becomes infinite. For that reasons, the 
construction of the classical solution along this branch is numerically difficult and we 
we had to stop it at A 3 « 0.4248. 

Nevertheless, it seems that the profile a(r) of the solution on the second branch 
is such that 

lim A3 ^ A . a(r) =| a BPST (r) | (46) 

a bp st being the profile of the charge- 1 instanton ( Jfl] ) for an appropriate value of 
the scaling constant k. The numerical difficulties prevented the evaluation of A3 but, 
according to Fig. 2, one can expect A3 ~ 0.42. 

The solutions were constructed with the subroutine COLSYS || (see Appendix 
of |I(J for a short description) and with a high degree of accuracy: typically with an 
error less than 10 -8 . 

To finish this Subsection we mention that the pattern of solutions presented above 
for A2 = seems to persist for A2 > 0. For instance, for A2 = 1 we find A3 ~ 0.14, 
i.e. a much lower value than in the case A2 = 0. More details are given in the next 
Subsection. 

4.2 Solutions with a(l) = a(oo) = 1 with A 2 > 0, A 3 > 

In this Subsection we present numerical results for solutions with the same asymp- 
totics as in Subsection |4.1| above, but (a) holding A2 fixed at A2 = 0.1 and varying 
A 3 , and (b) holding A 3 fixed at A 3 = 0.1 and varying A 2 . 

As in Subsection [O] above, the solutions do not appear to persist for arbitrarily 
large A3 (when A2 = 0.1), and arbitrarily large A2 (when A3 = 0.1). Unlike in 
Subsection [O] however, we have not endeavoured to find accurate critical values for 
the A2 and A3, respectively. The general features of the solutions remain unchanged 
whether or not A2 vanishes. 

The action versus A3 (with A2 = 0.1) is plotted on Fig. 3. The action versus A2 
(with A3 = 0.1) is plotted on Fig. 4. In both cases we note the remarkable feature 
that the action rises to just above the BPST instanton action |. This feature is shared 
with the restricted model with Ai = A3 = model studied in Ref. [|TJ. 

Fig. 5 is the analogue of Fig. 2, where the progress of r min and a m i n , the position 
of the minimum and the value of the minimum of the function a(r), are plotted against 
increasing A2 (with A3 = 0.1 fixed). Figs. 2 (resp. Fig. 5) describes the manner in 
which the solutions disappear as the value of the A3 (resp. A2) approaches a critical 
value. (The analysis confirming the existence of two distinct branches is given only 
in Fig. 2.) 

In addition to the above results, we have made a study of the generic system (|l|) 
with the potential (|j) included. It turns out that decoupling the potential (|j) results 
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Figure 4: A2 dependence of the action for A2 = 0.1. 
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Figure 5: A2 dependence of r m and a m for A3 = 0.1. 
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in no appreciable qualitative changes in the values of the action, for all values of the 
coupling constants A2 and A3. This property of the present model is shared by the 
usual (ungauged) Skyrme model 0. 

4.3 Solutions with a(0) = — a(oo) = 1 

Should solutions of this type exist, their Chern-Pontryagin charge would be nonzero. 
As we shall see below, the only such solutions are those with trivial S 4 valued field, 
i.e. only the pure YM Q instantons. 
The set of boundary conditions 

a(0) = l , a(oo) = -l , /(0) = tt , f(oo) = , (47) 

pertain to unit Chern-Pontryagin charge for the spherically symmetric configuration 
at hand, and provide a natural alternative to the solution discussed above. 

Although we could numerically exhibit regular configurations obeying these con- 
ditions and with an action slightly higher than 4/3, we have not succeeded in con- 
structing a solution of the equations with these conditions. We used both COLSYS 
and a relaxation method to find solutions satisfying the boundary conditions ([47|) 
and all the numerical results accumulated leads us to the formulation of the following 
statement: the solutions of the equations of motion obeying (f|7|) are constituted by 
the functions 

1,2 2 

< r ) = ~gr^ > keU , /(0) = , f(r)=n ,r>0 (48) 

considered in the limit k — > 0. The limiting configuration, which is determined in 
terms of step functions, has action 5 = 4/3. 

We now give an analytic construction to show that no such smooth solutions 
exist. To do this we will consider a suitably chosen one parameter family of field 
configurations and show that in the limit where the parameter vanishes, the field / 
becomes a step function shrunk near the origin and that the action reduces to the 
action of the instanton. 

First of all let us write the expression for the action after performing the scaling 
r = ax 

a 2 sin 2 / ( 2 q 2 sin 2 / ^ A 3 ( a 2 sin 2 A 2 ( 2 a 2 sm 2 f \ \ 3 

We then consider the following configuration: for a we take the usual instanton 
solution and we call xq the point where a(xo) = 0. For / we take 

/ = 7r if x < rr - I (50) 
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/ = (zo + |-aO^ if x - | < a; < x + | (51) 

/ = if x>x + |. (52) 

We then notice that the support of the action density for the first 3 terms is the 
interval [xq — |, xq + |]. Moreover in that interval we can write 

a = K{x — xq) + 0((x — xo) 2 ) (53) 

where K is a constant. Defining 

Sa = C-Y + i^ 1 ) 2 ) rHr (54) 



we can write 



1 r x +e/2 ( 

S = S a + - 

2 Jx -e/2 { 



a 2 X,x 3 (^ + - 2 K\x-x fsm 2 f } 



( K 2 it 2 K 2 \ 

+ A 2 x 3 \ —^{x- xq) 2 sin 2 /(— + —(x- x ) 2 sin 2 /) 
y x A e z x z J 

+ ^fx 3 (^-{x - x )) ^7r 2 + ^(x-x ) 2 sin 2 /) jrfx (55) 

We can the replace sin / by 1 and perform each of the integrals to the lowest order 
in e which gives 

S < S a + a 2 ^x n 2 + K 2 x e 3 ) + \ 1 a 2 O(e) + ^K 2 (^ 2 e + — ) 

+ «^ 3 > + ^ a ' 4 ((^W e3 + trw 7 ) + > (f5) (56) 

choosing <x = e we have 

lim S = S a (57) 

and in that limit the field / becomes singular, showing that there are no regular 
solutions with this boundary conditions. This is indeed what we observed when we 
tried to compute such solution numerically. 



5 Summary and discussion 

The coupling of non-linear sigma models to gauge fields often leads to sets of classical 
equations whose solutions obey various types of critical phenomena like bifurcations 
and/or pairs of solutions terminating into a cusp. The classical equations associated 
with the Lagrangian ([!]), in the spherically symmetric Ansatz, are of this type. These 
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solutions seem to follow the same pattern, irrespectively of the different Skyrme terms 
added, i.e. these patterns seem to be independent of the dynamical details. 

In this paper, we have studied the classical solutions ensuing from the Lagrangian ([!]) 
for three different sets of the two independent coupling constants (A2, A3). Inspite of 
the fact that our analysis in Section || (specifically in Subsection |2~T| ) leads to the 
establishment of a topological lower bound on the action provided that both A2 and 
A3 be positive, we have found that in fact solutions persist at A2 = 0. A similar sit- 
uation arises in the three dimensional Skyrme model augmented by a sextic Skyrme 
term. In that case too, when the usual (quartic) Skyrme term is decoupled, thus 
invalidating the topological lower bound, the solution persists [XT] notwithstanding. 
This most probably means that our (Bogomol'nyi type) analysis in Section Q is not 
far reaching enough for the model at hand. For example, in the case of Hopf solitons, 
there exists no Bogomol'nyi type lower bound on the energy, but instead one finds 
that a bound nevertheless can be established in terms of Sobolev type inequalities [|12[ . 



We have not carried out the appropriate analysis here, but expect that this can be 
done. Accordingly we have treated the simplified model (H) with A2 = as legitimate 
and have carried out the detailed analysis of exhibiting the cusp structure alluded to 



in the previous paragraph, for that model in Subsection |4.1| , which we summarise in 
the next paragraphs. 

For the A2 = model, clearly the two branches of solutions (the ones with non 
trivial S A valued matter field) terminate into a cusp at A3 = A3. This is a typical 
situation met in catastrophy theory. 

The spherically symmetric (BPST) instanton of the pure Yang-Mills theory [|J 
plays a major role and behaves as an attractor (at least when one of the coupling 
constants approaches a certain value) of the solution which excites both the matter 
and the gauge fields. It is very likely that the second branch of solutions bifurcates 
from the BPST-branch at the critical value A3 = A3 < A3. However, due to the 
absolute value in the limit P^|), the bifurcating solution does not occur in a standard, 
i.e. continuous way. 

The qualitative features of the instanton of the A 2 = model just described were 
confirmed also in the generic model with non- vanishing A2 and A3, in Subsections \L.2[ 
where the cusp resulting from the existence of two distinct branches was not explicitly 
constructed. 

In Subsection |4.3| , we verified that there existed no instantons with unit Pontryagin 
charge in this model, irrespective of the value of S 4 — *> S 4 winding number. This is 
important since it tells us that the zero Pontryagin charge instantons of this model 
are not the analogues of the triangle anomaly, and hence that the nonperturbative 
quasiclassical effects they describe must be given a new physical interpretation. (We 
return to this in the last paragraph.) 

Before alluding to the possible physical relevance of the model, we note that 
non-vanishing Pontryagin charge instantons can readily be constructed by changing 
the model to feature a symmetry breaking potential (|3|) as opposed to (f|). We 
have presented the simplest such model in Subsection |2~2] , but did not carry out a 
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numerical study in that case because the model involved was rather remote from 
known physically relevant models. 

In short, we have seen that the system consisting of the interacting YM and 
0(5) sigma models supports instantons with vanishing Pontryagin charge, which do 
not describe quasiclassical effects analogous to the triangle anomaly, but which have 
the S* 4 — > S* 4 winding number as the topological charge. Besides, the gauge group 
for this model is not that of the Standard Model. On the other hand, it is quite 
straightforward to construct an 0(5) model interacting with the 50(3) x SO (2) YM 
system that supports such instantons, by adapting the analysis of Subsection § to 
that case. (That remains a future project.) Moreover, the number of independent S* 4 
valued fields is equal to four, just as the number of the Higgs doublet in the Standard 
Model, thus an £0(3) x SO (2) gauged 0(5) model added to the Standard Model, 
could be regarded as a complicated low energy version of the latter, whose (axially 
symmetric) instantons can describe new nonperturbative effects. In this sense, the 
SO (4) gauged 0(5) model studied here can be regarded as a prototype of a physically 
more relevant model. 

Acknowledgements: This work was carried out in the framework of the TMR 
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